Large-mode-number magnetohydrodynamic instability driven by sheared flows 
in a tokamak plasma with reversed central shear 
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The effect of a narrow sub- Alfvenic shear flow layer near the minimum g m i n of the tokamak safety 
factor profile in a configuration with reversed central shear is analyzed. Sufficiently strong velocity 
shear gives rise to a broad spectrum of fast growing Kelvin-Helmholtz (KH)-like ideal magnetohy- 
drodynamic (MHD) modes with dominant mode numbers m, n ~ 10. Nonlinear simulations with 
finite resistivity show magnetic reconnection near ripples caused by KH-like vortices, the formation 
of turbulent structures, and a flattening of the flow profile. The KH modes are compared to double 
tearing modes (DTM) which dominate at lower shearing rates. The possible application of these 
results in tokamaks with internal transport barrier is discussed. 



Shear flows can give rise to the well-known Kelvin- 
Helmholtz instability (KHI) which, in magnetohydrody- 
namic (MHD) systems, has to compete with field line 
tension, magnetic pressure and magnetic sheari Re- 
gions in the plasma which are particularly susceptible 
to KHI are therefore those where it can couple with 
MHD instabilities. A system where such coupling is ex- 
pected to play an important role, and which motivated 
the present study, is the reversed-shear (RS) advanced 
tokamak configuration, which is a strong candidate for 
achieving self-sustained nuclear fusion conditions. It has 
a non-monotonic toroidal current profile [Fig[TJa)] and 
the magnetic shear, s = (r/q)dq/dr, changes its sign at 
some radius r m j n . Here r is the radial coordinate (mi- 
nor radius of the torus), and q(r) is the safety factor 
profile measuring the field line pitch [FigQJa)]. In the 
region with negative shear (r < r m ; n ) pressure-driven 
ideal MHD ballooning modes are stable, so that higher 
pressures can be achieved in the core. Strongly sheared 
poloidal flows may further improve confinement in RS 
tokamaks. There is evidence that such zonal flows play an 
important role in the formation of an internal transport 
barrier (ITB) near r m i n ££- MHD instabilities are often 
observed in ITB discharges.^ While there is evidence 
that certain ITB equilibria with strong pressure gradi- 
ent in the region of weakly positive shear may be stable 
to ideal MHD ballooning modes low shear around q m [ n 
may allow pressure-driven infernal modes to developi^ 
Recently, non-curvature-driven instabilities have also at- 
tracted attention!*^ The destabilization of an ideal elec- 
trostatic KHI in a region with low magnetic shear was 
demonstrated in Ref. [0]. KHI coupled to parallel mag- 
netic perturbations for s = was considered in Ref. Il2l . 

The scenario considered in this Letter consists of 
a radially localized shear flow layer near q m i n in an 
equilibrium which is also unstable to a broad spec- 
trum of current-driven resistive double tearing modes 
(DTM) ; 13 ' This may occur whenever q m - m passes 
through a low-rational value4£ Shear flow is known to 
stabilize DTMs through its tendency to decouple the 
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FIG. 1: Equilibrium profiles of (a) the safety factor q and 
current density j z o, and (b) the angular frequency uj-bq and 
vorticity u z o. q s = 2 resonances are indicated by dotted lines. 

gradient may destabilize DTMs . 17 ' 18 Thus, a dynamic in- 
teraction between DTMs and shear flows is expected^ 

In this Letter, we describe in cylindrical geometry how 
a broad spectrum of linear DTMs is stabilized by in- 
creasing velocity shear and show that, above a certain 
threshold, fast growing electromagnetic KH-like modes 
with dominant high mode numbers m ~ 10 appear. The 
linear instability of both DTMs and KH-like modes is 
analyzed and first nonlinear simulation results are pre- 
sented for a case with finite resistivity. The latter shows 
magnetic reconnection, the formation of turbulent struc- 
tures, and a flattening of the flow profile. 

For simplicity, we use the well-known reduced magne- 
tohydrodynamic (RMHD) model in cylindrical geometry 
(r, d, z).— In the limit of zero pressure, the RMHD equa- 
tions for the magnetic flux ip and the vorticity u z are 



dtif) + v • V?/> = —vjz 
(d t u z +vVu z ) = B Vj z 



(1) 
(2) 



Here, p m is the mass density, 77 the resistivity and v the 
viscosity. Current density j z and vorticity u z are related 



resonances ; ' In the nonlinear regime, locking between to ip and <fi through fioj z 



-Viip and u z 



finite-size magnetic islands and flattening of the velocity spectivcly. In the following, time t is measured in units 
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FIG. 2: (a) Linear growth rate 7n n and (b) rotation frequency ui^/m for the modes (m,n) = (2,1), (6,3) and (12,6), in 
dependence of the core rotation parameter G. The effective and specific shearing rates, A/ m and A/, are plotted as dash- 
dotted lines. In (c) spectra of linear growth rates 7n n (m) are shown for three cases with d = 0.02: G = stationary, G = 0.75 
stabilized, G = 1.5 destabilized. A spectrum obtained with a narrower shear layer is also shown (d = 0.01, G = 0.8). 



of the poloidal Alfven time th p = \/MoPma/i3o (with 
B being the strong axial magnetic field) and the radial 
coordinate r is normalized by the minor radius of the 
plasma, a. Resistive dissipation is measured by the mag- 
netic Reynolds number S'fjp = t^/thp (with t v = a 2 /j,o/n 
being the resistive diffusion time) and viscous damping 
by the kinematic Reynolds number i?en P = ci 2 /vtyi v . 
The perturbed fields are Fourier transformed as ft. = 
J2 m n ^m,«( r ) cxp(imi} — inz / R ) (with R being the ma- 
jor radius of the tokamak), where m/n = q s = 2 (single- 
helicity perturbation). The model and numerical meth- 
ods are described in detail in Refs. [3 and HO- 

The equilibrium configuration used is shown in Fig. [T] 
The safety factor has two q s = 2 resonant surfaces located 
a small distance D12 — r S 2 — r s \ = 0.066 apart (r sl — 
0.381, r S 2 = 0.447). The magnetic shear at these surfaces 
is = —0.10 and S2 = 0.12. The poloidal velocity profile 
v$o has a shear layer near r m j n and is given by 

= v# /r = fi G [C\ - tan -1 ((r - r )/d)] /C , (3) 

with Ci = tan _1 [(l - r )/d] and C = C± - 
tan _1 [(— ro)/d]. The standard case for uj$ is shown in 
Fig. QJb) . The shear layer is located near the inner reso- 
nance (ro = 0.4) and has a width (2d = 0.04) compara- 
ble to D12. The parameter G controls core rotation and 
thus the shearing rate between the resonances. We define 
the specific shearing rate as A/ = (cjq2 — ^oi)/27r [with 
u>0i = Lo$o(r S i)]. The effective shearing rate for the mode 
(m, n) is Af m = mAf. Since the equilibrium flow is not 
included self-consistently, only sub-Alfvcnic velocities are 
considered (Slo = 0.01, \G\ < 2). The characteristics of 
the shear flow are similar to those reported in Ref. 0. 
Both the static MHD equilibrium and the shear flow are 
separately unstable: the former to resistive DTMs^ and 
the latter to the ideal KHI>i the Rayleigh criterion for 
KHI is fulfilled [u' z0 (ro) = 0] and the existence of an in- 
stability threshold ™ ma x, is demonstrated shortly. 



For the shear-now-driven modes described in this Let- 
ter the vorticity gradient term p m (im / r)v r u' zQ (u' zQ — 
ditzo/dr) in Eq. ^fy is crucial. This term is known to be 
the origin of KHI.— In accordance with the criterion given 
by Eq. (23) in Ref. [l(| the presence of resonant surfaces 
is found to be crucial for instability. We find that low-m 
modes even require two nearby resonances and that mul- 
tiple branches of unstable modes (ideal and resistive) can 
be destabilized simultaneously. 

Results obtained by solving the linearized RMHD 
equations as an initial-value problem are presented in 
Figs. [2] and [3] In Fig.[2ja) the linear growth rates 7r m and 
in Fig. EJb) the linear frequencies wu n of several modes 
are plotted as functions of G. Spectra of unstable modes 
are shown in Fig.[5tc) for several values of G. The linear 
mode structures are shown in Fig. [3] The data were ob- 
tained with Sh p = f 7 and i?en P = f 11 , so the viscosity 
effect is negligible^ 

The 7iin(G) curves in Fig. [5fa) are symmetric around 
G = (no flow), so only G > is considered. Since the 
distance D\2 is small, the resonances couple strongly and 
a broad spectrum of resistive DTMs with dominant high 
poloidal and toroidal mode numbers m and n are found 
for G = (cf. Refs. [Hil. In the range < G < 1.0, the 
effect of the shear flow on DTMs is primarily a stabilizing 
one. However, it can be seen in Fig. Ela) that this stabi- 
lization may stagnate and even reverse. For instance, the 
growth rate of (m, n) = (12, 6) mode rises briefly around 
G ~ 0.2 before descending further. The origin of this 
behavior is not yet understood. The growth rates reach 
a minimum around G m - m = 1.0 and increase for G > 1.0 
almost linearly with G. When a smaller velocity gradi- 
ent (oc G/d) is used, or when the inter-resonance distance 
D12 is increased, a larger value is obtained for G m ; n (and 
vice versa). Without resonant surfaces ((7 m i n > 2) no 
unstable modes with helicity m/n = 2 arc found. 

The frequencies wr m (G) in Fig. [DJb) undergo a shift 
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around G m G mul . For G < G m - m the rotation at the 
outer resonance determines the frequency of the mode. 
This is typical for DTMs (since typically si < S2) and 
indicates that the current drive dominates in this regime. 
For G > G m in the frequency equals that at the inner res- 
onance, which is closer to the shear layer. Hence, G = 
G m in is identified as the transition point between pre- 
dominantly current-driven and shear- flow-driven modes. 
With a shear layer closer to the outer resonance there is 
no frequency shift. 

A comparison between the growth rates and the ef- 
fective shearing rates in Fig. [2](a) shows that 7u n (G) in- 
creases with G at least as fast as Af m . This suggests 
that some of the coupling between the resonances, which 
is important for DTMs, may pertain even for G > G m i n . 
For the (2,1) mode we have observed this effect by com- 
paring the results of the RS configuration with results 
obtained with a monotonic q profile with a single q s = 2 
resonant surface and otherwise similar conditions: in the 
former case (RS) the shear-flow-driven (2, 1) mode is un- 
stable [cf. Fig. EJ^a)], in the latter case (monotonic q) it 
is stable. In contrast, higher-m modes such as (6,3) and 
(12,6) have similar growth rates in both cases, so for 
these the coupling between the resonances does not seem 
to play a crucial role. 

The full spectra of linear growth rates 7n n (m) are 
shown in Fig. EJc). For the same shear layer width as 
in Figs. [2]^ a) and (b) (half width d = 0.02), the spectra 
are plotted for three values of G: no shear flow (G = 0), 
stabilized (G = 0.75), and destabilized (G = 1.5). The 
dominance of modes with m > 2 modes can be clearly 
seen. The DTM spectrum (G — 0) has m pea k = 6 as 
its fastest growing mode and the last unstable mode is 
"T-max = 18. In the stabilized case (G = 0.75) m pea k 
has decreased to 4. In addition, a second branch of 
weakly unstable modes is present, with m pea k = 18- For 
G = 0.75 the frequencies of modes on both branches are 
determined by the rotation of the outer resonance. In the 
destabilized case (G = 1.5), the modes on the dominant 
branch around m poa k = 10 grow rapidly and their fre- 
quencies are up-shifted, implying rotation with the inner 
resonance. The growth rates on the secondary branch are 
almost the same as for G = 0.75 and the frequencies are 
still determined by the rotation of the outer resonance. 

In Fig. [2]Jc) a case with a smaller shear layer width 
d = 0.01 is shown as well. Compared to the d = 0.02 
case, here the dcstabilization of shcar-flow-driven modes 
has a lower threshold G m ; n and the spectrum is broader: 
mpeak = 14 and m max = 26. Indeed, it is found that 
ra max is controlled by the parameter d. while G seems 
to be irrelevant. This is similar to KHI which possesses 
an upper threshold for the wave number k$ = m/r of 
unstable modes:— for 7 > it is required that k$d < 0(1). 
Here, for d = 0.02 and tq = 0.4 the dominant branch has 
Wmax = 18 [Fig. |2£c)], so that k$d = 0.9 < 1. Modes on 
the secondary branch do not seem to be subject to the 
same threshold. 

Calculations with magnetic Reynolds numbers in the 
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FIG. 3: Linear eigenmode structures for (m, n) = (10, 5) [(a) 
and (b)] and (20, 10) [(c) and (d)]. The flow parameters are 
to = 0.4, d = 0.02 and G = 1.50. Flux ip and potential <j> 
are decomposed into magnitude and complex phases as tp = 
j^e 10 ^. Vertical dotted lines indicate resonant surfaces. 



range 10 < Sh p < 10 15 reveal that shear-flow-driven 
modes (G > G mm ) on the dominant branch grow inde- 
pendently of Sn p . The ideal character of the instability 
is reflected in the linear eigenmode structures shown in 
Fig. Ela) and (b) for (m,n) = (10,5): the flux perturba- 
tion switches signs at the resonances, as can be inferred 
from |^|(r S j) = and the phase shifts Aa(r S i) = ±7r 
[Fig.|3Ja)]. The potential <f> [Fig.[3{b)] is non-zero at the 
resonances and its complex phase exhibits a jump at the 
inflection point of the velocity profile, where u' z0 = 0. 
Resistivity independent modes are also observed with a 
single resonant surface and otherwise similar conditions. 

The typical structure of modes on the secondary 
branch is shown in Fig.[3]Jc) and (d) for (m, n) = (20, 10). 
The perturbation is localized between the shear layer (tq) 

and the outer resonance (r - ^). At r — r S 2, the flux ip is 
non-zero and the potential <j) has a discontinuity, sug- 
gesting a tearing-type instability. This mode grows most 
rapidly for 5h p ~ 10 7 , while the growth rate is reduced 
for larger and lower resistivity. 

Further calculations show that the growth rates on the 
dominant branch are rather insensitive to the exact lo- 
cation of the shear layer between the resonances. This 
indicates that shear-flow-driven modes have a non-local 
character, in the sense that it is not sufficient to consider 
values of u>$o at the resonances only (cf. Ref . [Toh . 

Finally, nonlinear simulation results are presented in 
Fig. [4] (parameters settings: G = 1.5, d — 0.02, 5h p = 
10 7 , i?e Hp — 5 x 10 8 , single helicity q s = 2, 128 modes, 
dealiased). At the time of the first snapshot, Fig. 0|a), 
the nonlinear mode coupling has reduced the velocity gra- 
dient around tq = 0.4, while the profile in the surround- 
ings is steepened. The helical flux contours show islands 
around the inner resonance, with a dominant m = 2 
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FIG. 4: Snapshots of nonlinear dynamics at (a) t = 1200 and 
(b) t = 1770. Top: lj-b profiles. Dotted lines indicate the 
original resonant radii. Middle: q profiles. Bottom: helical 
magnetic flux contours [ip, = ip + r 2 /(2q B )] in the d-r plane. 

topology and high-m ripples. In the X-point regions, 
poloidally localized KHI-like vortices (only one shown 
due to 7r-periodicity) can be observed. Subsequently, the 
m = 2 islands undergo fragmentation and the structures 
become increasingly turbulent, as can be seen in the sec- 
ond snapshot in Fig.[4jb). 

Note that the shear-flow-induced perturbations ex- 
tend over the entire inter-resonance region (and beyond). 
Thus, one may expect a coupling similar to that in 
DTMs. Indeed, Fig.|4jb) shows island formation on both 
resonances while the central flux surfaces are still intact. 
Furthermore, in Fig. Q|b) a flattening of the flow profile 
in a wide region can be observed, while g m i n is still well 
below q s = 2. This shows that when resistivity is low, 



the primary saturation mechanism in the early nonlinear 
stage is likely to be the flattening of the flow profile. 

The observed magnetic reconnection is due to finite 
plasma resistivity rather than turbulence effects. How- 
ever, there are two instability driving mechanisms, the 
current gradient and the velocity gradient, and it is not 
clear how much each contributes at a given time. The 
observed flattening of the flow profile suggests that the 
current-gradient-driven instability (i.e., DTM) may play 
a more significant role in later stages. The study of the 
long-time evolution of this system requires turbulence 
simulations which go beyond the scope of the physical 
model and numerical code used here. 

In summary, we have analyzed the effect of sheared 
flows in a RS tokamak configuration when q m i n is below 
but close to a low-order rational value q s = m/n and a 
pair of nearby resonant surfaces has formed. The growth 
rates of DTMs are reduced, although the spectrum re- 
mains peaked at m > 2 even with strong shear flow. For 
sufficiently strong velocity shear DTMs are replaced by a 
broad spectrum of shear-flow driven KH-like modes with 
dominant high m ~ 10. In the nonlinear regime and 
with finite resistivity, magnetic reconnection occurs near 
ripples caused by KH-like vortices. The tendency to de- 
velop turbulent structures in the inter-resonance region 
and flattening of the velocity profile is observed. Further 
calculations with experimental data as input are required 
in order to identify the role of KH-like modes in ITB dis- 
charges. Note also that since high-m modes are involved 
the single-fluid MHD picture may be inaccurate. Our re- 
sults motivate further investigations with more realistic 
models to understand the properties and effects of KH- 
like modes in tokamak plasma with ITB. Since high-m 
modes are involved, we expect that similar instabilities 
may be found for other values of q s such as 3/2 or 7/3*^ 
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